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> We rigorously estabiish the vaiiditv of the equation, describing the evointion of one- 

^ j dimensional long wavelength modulations of counterpropagating wavetrains for a hyperbolic 



model equation, namely the sine-Gordon equation. We consider both periodic amplitude func- 
tions and localized wavepackets. For the localized case, the wavetrains are completely de- 

(N 

coupled at leading order, while in the periodic case the amplitude equations take the form 
O 

of mean-held (nonlocal) Schrodinger equations rather than locally coupled partial differential 

equations. The origin of this weakened coupling is traced to a hidden translation symmetry 



in the linear problem, which is related to the existence of a characteristic frame traveling at 



the group velocity of each wavetrain. It is proved that solutions to the amplitude equations 

^— > 

C$ dominate the dynamics of the governing equations on asymptotically long time scales. While 

the details of the discussion are restricted to the class of model equations having a leading 
cubic nonlinearity, the results strongly indicate that mean-held evolution equations are generic 
for bimodal disturbances in dispersive systems with 0(1) group velocity. 
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1 Introduction 

There has been considerable progress recently in demonstrating the validity of amplitude equa- 
tions as approximations to the governing equations of dissipative systems [2, 4, 12]. For the 
most part, this work has been restricted to systems for which the initial bifurcation is either to 
traveling waves or steady state solutions. We concentrate on extending the results of Kirrmann, 
Schnieder and Mielke [4] (KSM, hereafter) to the case of counterpropagating wavetrains, of 
which traveling and standing waves are special cases. KSM are among the few that have 
considered hyperbolic systems (cf. [3] also) in addition to the more common dissipative case. 
We follow their derivation and prove that solutions to the amplitude equations derived below 
remain close to solutions of the governing equation over asymptotically long time scales. As a 
model we choose the governing equation to be the one- dimensional sine-Gordon (SG, hereafter) 
equation. However, the structure of the proof is unchanged for any equation having translation 
and reflection symmetries in both space and time as well as a leading cubic nonlinearity. The 
translation symmetries imply that the imposition of periodic boundary conditions results in a 
system which is 0(2)-symmetric. When this property is combined with the fact that the group 
velocity is of 0(1), asymptotic scaling arguments indicate that the correct modulation equa- 
tions are of the mean-field type [5, 6, 7, 9, 10]. However, in related systems other arguments 
have lead to the conclusion that the coupling should be local [8]. In an effort to resolve this 
question, we first derive the mean-field nonlinear Schrodinger equations (MFNLS, hereafter) as 
the asymptotically exact amplitude equations for the SG equation. We then rigorously prove 
that the approximate solution represented by the amplitudes remains asymptotically close to 
a full solution to the governing equation on time scales of 0(e~ 2 ). Here e is a small parameter 
measuring the size of the disturbance. In addition, we show that similar estimates hold both 
for solutions in if 1 (i?), and also in spaces of periodic functions. We conclude with a discussion 
of how our results are relevant to studies of wave propagation in physical systems. 
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2 The Amplitude Equations 

We look for solutions u(x,t) : R X [0, oo) — ► R to the SG equation, 

d 2 u = d 2 u - sin-w . (1) 

If we restrict attention to solutions of the linearized problem with amplitude of 0(e) and 
periodic in x with period ^, we have: 

u L = e A+e i(kx - wt) + e A~ e i(kx+UJt) + cc , (2) 

where the A ± are 0(1) constants and cc denotes the conjugate of all preceding terms. The 
frequency uj is related to the wavenumber k by the dispersion relation 

cj 2 = k 2 + 1 . (3) 

u L is a valid solution to (1) in the limit that e —^ 0. Note that (1) has a translational 
symmetry in both x and t. Each wavetrain in u L breaks these symmetries, but each preserves 
a spatio-temporal symmetry in which translations in space are equivalent to translations in 
time. Together they break the remaining symmetry because the spatial translations are in 
opposite directions for a given temporal displacement, but the wavetrains are decoupled in the 
linear problem and so each retains its symmetry properties independent of the other. 

We continue to examine the linear case, but we now relax the periodic boundary conditions, 
allowing for a continuous wavenumber spectrum. The spatio-temporal translation symmetry 
of each wavetrain is broken by modulation of the envelopes in space and time. In this case, it 
appears that all of the symmetries of the problem are now broken. However, modulations of 
the wave envelope with infinitesimal wavenumber propagate at the group velocity, so that each 
envelope is time independent in the characteristic frame moving at its group velocity, c g = — 
for the SG equation. Consequently, there is a hidden spatio-temporal translation symmetry 
associated with each wave envelope, for which translations in space correspond to translations 
in time. It is important to note that this hidden symmetry is a property of the linear problem, 
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but it is a symmetry of the amplitudes and does not correspond to a symmetry of the full 
modulated linear solution. The net result of this argument is to establish that if we wish to 
describe modulated linear waves, then the amplitudes should be functions of the characteristic 
variables A ± = A ± (x ^ c g i). 

For finite but small e, the nonlinear terms introduce perturbations to the linear problem. 
This causes the amplitudes to vary in time and space, and it also causes the wavetrains to 
interact. However, since the perturbed solution must converge to the linear solution in the 
limit that e — ► 0, the amplitudes may vary only on scales which are asymptotically long in the 
limit, i.e., A ± = A ± (et, ex, e 2 t, ■ ■ •), and interactions between the wavetrains must take place 
at higher order in e. Due to the hidden translation symmetry, the amplitudes cannot have 
arbitrary dependence on the slow variables ex and et but must instead depend only on the 
slow characteristic variables x ± — e ( x T c g~t) a t leading order. Moreover, terms which break 
this hidden symmetry can only enter at higher order in e. Note that we can only make this 
argument because the hidden symmetry is a property of the linear problem. The approximate 
solution which is expected to be valid for some finite range of e should then have the form: 

u A = e(A + (x + ,r) + eB+ ( X + , X~ , r))e^-^ + e( A~ ( x ~ , r ) + eB~ ( X + , X~ , r))e^ + ^ 

+ cc , (4) 

where r = e 2 t is a super-slow time scale. The B ± are the leading order terms which break the 
hidden symmetry, and, like the A ± , they represent long wavelength modulations of the linear 
solution. As noted above, they must enter the expansion at 0(e 2 ) at the earliest because only 
nonlinear effects can break the hidden symmetry. There are no second harmonics because the 
nonlinearity enters at 0(e 3 ). We have not included a super-slow spatial scale, £ = e 2 x, because 
our results prove that it is not necessary in order to capture the leading order behavior for 
time scales of 0(e~ 2 ). 

Following KSM, we introduce an "improved" approximation, from which we will derive the 
amplitude equations: 

v A = e(A+ + eB+ yC**""*) + e{A~ + eB~ ) e i(kx+ ^ 
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l f 3 / ( A+ 3 3i(kx-wt) , A-3 3i(k X +wt)\ 

J ( A+ 2 A- J(3kx-Lvt) _|_ a+ a- 2 i(3kx+wt)\ 

2(9k 2 -uj 2 + iy ' 

+ 2(k2 _ ^ 2 (A+'l^e^^- 3 ^) + A+A- V^ 3 "*))} + cc , (5) 

where the overbar denotes conjugation. The cubic terms are designed to cancel the nonres- 
onant terms resulting from the expansion of the nonlinearity in the governing equation. For 
the current problem, these higher harmonics are nonresonant for all k ^ 0, but for other prob- 
lems this is not always true and the resonant harmonics must explicitly be introduced to the 
expansion at lower order. We next compute the amount by which v A fails to be a solution by 
substituting it into the governing equation and calculating the "residuum": 

p(e, x,t) = d 2 v A - d' 2 x v A + sin v A 

= e 3 [ - 2iud T A+ - Aikd x -B+ + (c 2 - 1)^+A+ - -\A+\ 2 A+ - \A-\ 2 A+\e i{kx -^ 
+e 3 \2iL>d T A~ -4ikd x+ B- + (c 2 g - l)d 2 x -A~ - -\A~\ 2 A- - |^+|2^-] e ^— «) 
+ cc+0(e 4 ). (6) 

The quantities in the square brackets are essentially the amplitude equations, but the evolution 
of both the first and second order amplitudes cannot be resolved from the two equations. 
The correct way to proceed is to separate the terms which depend on the counterpropagating 
characteristic variables into separate equations. This is equivalent to interpreting the amplitude 
equations as inhomogeneous equations for the B"* 1 and applying the consequent solvability 
condition (cf. Section 4 for further discussion). 

a-)A+ (7) 

a+)A~ (8) 

(9) 
(10) 

The a ± are functions of the super-slow time r only, and their definitions depend on the function 



2iujd T A+ -- 


= (cl-l)d 2 x+ A+-q\A+\ 2 


-2icod T A- = 


= (cl-l)dl-A--( 1 -\A-\ 2 


4ikd x -B+ = 


= -(\A-\ 2 -a~)A + 


4ikd x+ B~ -- 


= -(]A+\ 2 -a+)A- 
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space in which we seek solutions. If solutions are sought in the form of localized wave packets, 
then 

^(^ = 0. (11) 

If the amplitudes are forced to be periodic of period P l5 the a ± are defined as 

o ± (r)^^ FlA^rtfds. (12) 

-* l Jo 

The terms B^ have now decoupled from the amplitudes, which must evolve according to the 
MFNLS equations, (7) and (8). The solutions of (9) and (10) for B ± may be written explicitly: 



B ± (x\x-,T) = ^A ± ( X ± ,T)j* T (\A*(s,T)\ 2 -a*(T))ds, 



:i3i 



where £ = — oo in the whole-line case and £ = in the periodic case. Note that the B^ remain 
of 0(1) regardless of whether the A ± are localized or periodic. Note also that the B^ are 
determined by (9) and (10) only up to the addition of an arbitrary function 7 ± (x ± ,r). In 
writing (13), we have set 7 ± = 0. A different choice of functions may be physically relevant 
in some cases, but this presents no additional complications provided that the 7 ± (-,r) are 
bounded in the appropriate norm. 

3 Estimates on the Validity of the Amplitude Equations 

While the preceeding arguments are convincing, we now show that they can, in fact, be made 
mathematically rigorous. We begin by discussing the case of localized wavetrains, since the 
discussion in this case parallels that of KSM. Let u(x,t) be the solution of the SG equation. 
As KSM point out, the natural space in which to consider this problem is (u,d t u) £ Y = 
H 1 {R) X L 2 (R). The global well-posedness of SG in this space is a standard result. Similarly, 
by the results of [1] (section 6), one has global well-posedness of (8) and (9) (with a ± = 0) 
in if 1 (i?). Finally, consider B"* 1 which solve (9) and (10) (and for which we have an explicit 
representation in (13)). From (13) one has immediately that 
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As an immediate consequence of this estimate we have: 
Lemma 1 

/oo ^ 

^(ex + ec g t, ex - ec g t, e 2 t)\ 2 dx < y^ll^O, ^)WlAW(-, e 2 t)\\ 4 L2 , (15) 

while 

H^c-.-^OlU- < ^II^C-^oiU-IIA^ni, . (16) 

Remark 1 By the Sobolev inequality, ||^4 ± |||oo < 2||A ± || i 2||A ± || ff i. 

Note that these estimates allow us to bound the L 2 and L°° norms of our approximate solution 

v A . 

Corollary 1 Let ||^4||l 2 = max (||^- + IU 2 5 ||^-~||l 2 ) o,nd similarly for \\A\\ H i . Then for e suffi- 
ciently small, there exists a constant c(k,u>) such that \\v A \\ H i < c (^ 5 w )e 1/ ' 2 ||^4||ff 1 - 



< II 4±|| 2 || 4=F||2 /i 7 \ 

- ^u, l II \\L II \\L V ) 



Proof: The estimates of the lemma allow us to bound 

A (ex ± c g t)B (ex + c g t, ex — c g t, e 2 t) dx 

Thus, taking into account the scaling of the spatial variable, we find 

IM-,r)||i 2 < c{k,u) (e\\A\\ 2 L2 + e 2 \\A\\ 4 L2 + e 3 \\A\\l 2 \\A\\ m + e 4 \\A\\ e H1 ) . (18) 

For e sufficiently small, this is bounded by c(A;,cj)e||A||^i. The estimate on H^llff 1 is similar 
and we don't repeat it. rj 

Note further that as a corollary of these estimates on A ± and B ± and the fact that they 
satisfy (7) - (10), we obtain the following estimates of the residuum (6): 

Corollary 2 Suppose that A^ are solutions of (7) and (8), for which d"d^±A^ G 

C([0,T ], L 2 (R)) for a + (3 < 2. Then there exists a constant c > 0, depending on k, uj, and 

the norm of A ± and its derivatives of order 2 or less, such that for all < t < T e~ 2 , 

IKM,-)IU~ < « 4 (19) 

\\p(e,t,-)\\ L , < ce 7 ' 2 (20) 

\W,t,-)\\ m < ce 7 ' 2 . (21) 
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Remark 2 If one assumes that the initial data A ± (x,0) of (7) and (8) is sufficiently smooth, 
it is straightforward to show that d"d ±A ± £ C([0,T ], L 2 (Rj) using ideas like those used by 
Cazenave ([1], section 5) to prove well-posedness in H 2 . 

We now estimate the error R(x,t) = e~ 3 / 2 (u(x,t) — v A (x,tj), where u(x,t) is a solution of 
the SG equation. Given the estimates on v A , we can follow KSM almost word for word through 
these estimates. Note first that 

d 2 R = d' 2 R- R + a(e, t)R + e" 3/2 iV(e, t, R) - e" 3/ 2 p(e, t) , (22) 

where a(e, t) = (1 - cos(w^)), and JV(e, t, R) = e 3 | sin(w^ + 6e 3/2 R)R 2 , for some \0\ < 1. 

Remark 3 If we had chosen to consider a more general nonlinear equation than the SG equa- 
tion, say d 2 u = d 2 u — u-\-g(u), the only change in this expression would be that a(e,t) = g'(v A ), 
and N(e, t, R) = e 3 \g"{v A + 9e 3 l 2 R)R 2 . 

As KSM note, if we rewrite (22) as a first order system in terms of the variables (R, S) = 
(R,dtR), then the linear evolution 



R\ _ ( S 
S ) ~ \dlR-R 



M c) = ( .: J (23) 



preserves the norm \\{R,S)\\ 2 Y = J™Jd x R) 2 + R 2 + S 2 dx. Since \\v A \\ L ~> < Ce, \a(e,t)\ < Ce 2 
for all t, while N(e,t,R) is a smooth function from Y to L 2 with ||iV||£2 < e 4 ||-R||jyi. Thus, we 
can rewrite (22) as the integral equation: 

( 5(0 ) = G{t) im) + l G{s - t) ia(^)R(s) + ,- 3 / 2 iV( e , t,R)- e~^p(e, s) ) ds ' < 24 > 

where G{t) is the semigroup associated with the linear equation (23). Using the fact that G 
preserves the Y norm, standard estimates imply that 

\\(R(t),S(t))\\ Y < \\(R(0), 5(0))|| y + 6- 3 / 2 sup \\p(e,s)\\ H1 t 

0<s<t 

{C^WiRis), S(s))\\ Y + C 2 e 5 ' 2 \\(R(s), S(s))\\ 2 Y } ds , (25) 
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where we have used the fact that if ||_R(s)||#i < C, then || sin(f^ + 0e 3 / 2 R)\\l<*> < C 2 e, for e 
sufficiently small. Next presume that we can choose our initial amplitudes A ± (x ± ,0) so that 
||w(a;,0) — v A (x, 0)||#i < Ce 3 / 2 , and ||9 t M(a;,0) — d t v A (x, 0)\\l 2 < Ce 3 / 2 . Thus, we can assume 
that ||(i2(0), 5'(0))||y < C\. Similarly, using our estimate on \\p(e, s)\\ H i, from Corollary 2, 
we see that e -3 / 2 sup 0<s<t ||p(e, s)||#ii < CT . Thus, we can apply Gronwall's inequality as in 
KSM to conclude: 

Theorem 1 Suppose that A^ are solutions of (7) and (8), for which d"d^±A^ G 

C([0,T ], L 2 (R)) for a + (3 < 2. Then there exists e and C greater than zero, such that if 

< e < e the solution of (24) satisfies \\(R(t), S(t))\\ Y < C , for all < t < r e" 2 . 

From this one immediately concludes: 

Corollary 3 Under the hypotheses of the previous theorem, there exist positive constants e 
and C such that if e < e , then 

\\u(-,t)-v A (-,t)\\ m <C e 3 l 2 , (26) 

for all < t < T e~ 2 . By Sobolev's inequality, one immediately has 

\\u(-,t)-v A (-,t)\\ L ~<C e 3/2 , (27) 

for the same time interval. 

We now turn to a consideration of the problem with periodic boundary conditions. We 
consider solutions u(x,t) which are periodic with period P e , where P c and ^ are commensurate. 
This implies (because of the scaling of the spatial variables) that the amplitude functions A ± 
and B ± are periodic with period eP e , and hence we choose P e = e _1 Pi and Pi is an 0(f) 
constant. Define a norm 

IK,i)lly P „(0 = P" 1 [\u 2 + (d x u) 2 + (d t u) 2 } dx . (28) 

We also define the norms ||w(-,i)||£2 ( e ) = P~ x f ' u 2 dx, and by analogy Pp er (e). 
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In like fashion we will define the -ff* norm for A ± by 



H^^r)!!^ =P 1 " 1 / 1 {|A ± (x ± ,r)| 2 +|a x± A ± ( X ± ,r)| 2 }d X ± , (29) 

per J Q 

and analogously for L 2 r . 

Remark 4 The well-posedness of (7) and (8) follow immediately if we note that cf ± {t) = 
P ± (-' r )lli=„ = \\ A± (-i°)\\l\ r are independent of t and define A ± = i±exp(±^). The A ± 
then satisfy uncoupled NLS equations. 

Now consider B^ defined by 

B ± (x\x-,T) = ^A ± ( X ± ,T)j o XT (\A*(^T)\ 2 -a*(T))dt. (30) 

Note that the periodicity of A ± and the presence of the counterterm a T ensure that B^ is 
periodic in both variables as required. Similarly, one can repeat the estimates leading up to 
Lemma 1 and one finds: 



Lemma 2 

/o 



f p ' 1 

P- 1 ^(ex + ec g t,ex - ec g t,e 2 t)\ 2 dx < — — - 1 1 ^± ( - , e 2 ^) | |i 2 (e) ||A ± (-, e 2 i 






(31) 

1 

4IA-I 

(32) 



|^(-,-,^)IU~ < ttuII^C-^OIU-II^C-^OIIImo 



One can now proceed to estimate our approximate solution v A exactly as we did in the 
whole-line case. We find 

Lemma 3 Let IIAIIh-i = max(||A + ||H-i , ||A _ ||h-i ), and similarly for II All roo. Then 

II ii per ^'1 ii per' II II per' 7 O J II | | u 

Pe' 1 I ' {\v A {x,t)\ 2 +\d x v A {x,t)\ 2 )dx < C(k,u){(e\\A\\ HL +e 2 \\A\\ Hi ) 2 



per s 



-^\\MU\Mn i .M + \\A\\ i Hiar + \\A\\U)}. 

(33) 
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Remark 5 There is a slight difference here between the periodic and whole-line problems which 
concerns the estimate of \\A\\ L ^ . In order to estimate this, choose x to be points such 
that I^O^)! 2 = a ± . Then |A ± (a;)| 2 - a ± = 2 f*± A ± (£)d x ±A ± (£) d£, so we conclude that 

II /t ± 112 < n -±_|_0D||/1||2 

/I Too ^ V \ 1 /T 1 

The important point is that although somewhat more complicated than than the corre- 
sponding estimate in the whole-line situation, the right hand side of this inequality is still 
bounded, independent of e for all < t < T . 

In like fashion we can derive estimates on the residuum (6). We find: 

Corollary 4 Suppose that A^ are solutions of (7) and (8), for which d^d^A^ £ 

C([0,T ], T 2 r ) for a + (3 < 2. Then there exists a constant c > 0, depending on k, uj, and the 

norm of A ± and its derivatives of order 2 or less, such that for all < t < T e~ 2 , 

IKM,-)IU~ < « 4 (34) 

IMM,-)||*i„ < C6 4 . (35) 

Once again, we emphasize that these estimates can be chosen independent of time, for < t < 
T e- 2 . 

We complete our analysis of the periodic case by estimating the error R per (x,t) = 
e~ 2 (u(x,t) — v A (x,t)). (Note the change in scaling with respect to the whole-line case. This 
results from the change in the way the norms scale with e in the periodic case.) 

where as before a(e,t) = (1 — cos(t^)) and N per (e,t,R per ) = e 4 |sin(f^ + e 2 9R per )R 2 r . As 
before we note that the semi-group G per associated with the linearized equations (36) (subject 
to periodic boundary conditions) preserves the norm ||(i? per , <5 , per )||y • Hence, as before we 
can estimate (36) by rewriting it as an integral equation and we find: 

\\(R per (t),S per (t))\\ Ypsr < \\(R per (0),S per (0))\\ Ypsr + e" 2 sup \\p(e,s,-)\\ H i (e) t 

0<s<t 
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t 

2 







{Ce 2 \\(R per (s), S per (s))\\ Yp „ + Ce*\\(R per (s), S per (s))\\ 2 Yp J ds . 

(37) 

Just as in the whole-line case, we choose our initial amplitudes so that 
||(-R per (0), 5 , per (0))||y er < C. Thus by another application of Gronwall's lemma, we conclude: 

Theorem 2 Suppose that A^ are solutions of (7) and (8), for which d^d^A 42 G 

C([0,T ], L 2 r ) for a + (3 < 2. Then there exists e and C greater than zero, such that if 

< e < e the solution of (37) satisfies \\(R(t), S(t))\\ Ypsr < C , for all < t < T e~ 2 . 

As before, we can immediately conclude from this that ||w(-,i) — v A (-,t)\\ H i ^ < C e 2 , for 
all < t < T e~ 2 . But then if we choose x so that (u(x ,t) — v A (x ,t)) 2 = P~ f '(u(x,t) — 
v A (x,t)) 2 dx < C e 4 , we have 

(u(x,t) - v A (x,t)) 2 = (u(x ,t) - v A (x ,t)) 2 

+2 / (d x u(C,t)-d x v A (C,t))(u(C,t)-v A (C,t))dC 

■)x 

< C e 4 + f{(d x u(^ t) - d x v A {£, t)) 2 + (u(i, t) - v A {£, t)) 2 } dti 

< C e 4 + 2P e \\u - v A \\ 2 H}er(e) < C e 4 + 2P 1 e 3 . (38) 

Thus, we have proven the following corollary: 

Corollary 5 Under the hypotheses of the preceeding theorem, there exist positive constants e , 
C\, such that if < e < e then the solution of (36) satisfies ||w(-,i) — ^(^OIU 00 — Cie 3 ^ 2 for 
allO<t< T e~ 2 . 

4 Discussion 

We begin by considering alternative approaches to eliminate the 0(e 3 ) terms from the residuum 
(6), and describe why we feel our method is the most appropriate. The quantities in the square 
brackets are recognizable as NLS equations with a cubic cross term and, effectively, a slow time 
derivative of the symmetry breaking B 41 terms. These derivatives may well be regarded as the 
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"problem" terms. It is tempting to argue that they might be neglected because they are higher 
order terms, either by breaking the hidden symmetry at leading order or simply setting the 
B ± to zero. In the former case, one might try to define A ± = A ± (et,ex, ■ ■ •), but the hidden 
symmetry is restored at 0(e 2 ) via the evolution equations 

d xT A ± = 0. (39) 

The right-hand side of these equations are identically zero due to the absence of a quadratic 
nonlinearity in the governing equation. If the alternative approach is taken and the B^ are 
set to zero, the amplitude equations take the form of locally coupled NLS equations (i.e., 
(7) and (8) with a ± = \A ± \ 2 ). Considering the evolution equation for A + , it is independent 
of the characteristic variable x~ except for a parametric variation of the term |A~| 2 A + , and 
similarly for the A~ equation. Grouping terms depending on x + an( i X~ on opposite sides 
of the equations, it is apparent that the quantities \A ± \ 2 depend only on r and not on the 
characteristic variables x ± ? at least as long as the amplitudes do not have zeroes. Consequently, 
the amplitudes take the form A^ = a ± (T)e 1 ^' ^ x ,T \ and the locally coupled NLS equations 
are found to exclude the possibility of amplitude modulation. While one must entertain the 
possibility that this is the physically relevant description, there are two compelling arguments 
against it. Mathematically, the process of setting the B ± to zero contains the unjustified 
assumption that the terms which break the hidden symmetry (and destroy the existence of 
the characteristic frames) appear at 0(e 3 ) and not at the first opportunity as one generically 
expects. Physically, the exclusion of amplitude modulation is contrary to the existence of the 
Benjamin- Feir instability, which is a well-documented feature of wavetrains in physical systems. 
The results presented above generalize those of KSM in several ways. Most importantly, we 
have addressed the case of counterpropagating wavetrains and the case of periodic boundary 
conditions. The details of the proof restrict our results to governing equations having a leading 
cubic nonlinearity, thereby excluding most systems of hydrodynamical interest. It is likely that 
these results will hold for systems having quadratic nonlinearities, provided that the possibility 
of resonant triads is excluded, but it should be noted that it is not a trivial exercise to prove this 
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rigorously [11]. Moreover, we have shown that the mean-field terms in the amplitude equations 
originate with a hidden symmetry of the linear problem. Consequently, it is unlikely that the 
structure of the nonlinearity will have a substantial effect on the nature of the mean-field 
coupling of the wavetrains. 

It has also been demonstrated recently that the asymptotically exact equations governing 
small-amplitude broad-band disturbances (the so-called AEZ equations) are similar to the 
well-known Zakharov equations [13], but the AEZ equations differ in that they contain a time 
average [9]. At the same time it was shown that the averaging is responsible for generating the 
mean-field terms in, for instance, the MFNLS equations derived here; the original Zakharov 
equations without the average yield locally coupled NLS equations. Consequently, we interpret 
the present results as strong, though not conclusive, evidence for the generic nature of mean- 
field coupling for both broad-band and nearly monochromatic small-amplitude disturbances. 
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